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Field Induced Multiple Reentrant Quantum Phase Transitions in Randomly 
Dimerized Antiferromagnetic 5*= 1/2 Heisenberg Chains 
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The multiple reentrant quantum phase transitions in the 5 = 1/2 antiferromagnetic Heisen- 
berg chains with random bond alternation in the magnetic field are investigated by the density 
matrix renormalization group method combined with the interchain mean field approximation. 
It is assumed that the odd-th bond is antiferromagnetic with strength J and even-th bond 
can take the values Js and Jw ( Js > J > Jw > 0) randomly with probability p and 1 — p, 
respectively. The pure version (p = and p = 1) of this model has a spin gap but exhibits 
a field induced antiferromagnetism in the presence of interchain coupling if Zeeman energy 
due to the magnetic field exceeds the spin gap. For < p < 1, the antiferromagnetism is 
induced by randomness at small field region where the ground state is disordered due to the 
spin gap in the pure case. At the same time, this model exhibits randomness induced plateaus 
at several values of magnetization. The antiferromagnetism is destroyed on the plateaus. As a 
consequence, we find a series of reentrant quantum phase transitions between the transverse 
antiferromagnetic phases and disordered plateau phases with the increase of the magnetic field 
for moderate strength of interchain coupling. Above the main plateaus, the magnetization curve 
consists of a series of small plateaus and the jumps between them. It is also found that the 
antiferromagnetism is induced by infinitesimal interchain coupling at the jumps between the 
small plateaus. We conclude that this antiferromagnetism is supported by the mixing of low 
lying excited states by the staggered interchain mean field even though the spin correlation 
function is short ranged in the ground state of each chain. 

KEYWORDS: random quantum spin chain, DIVIRG, disorder induced order, field induced order, randomness 
induced plateau, reentrant phase transition 



1. Introduction 

In the recent studies of one-dimensional quantum spin 
systems, the exotic quantum phases induced by the 
strong magnetic field have been attracting broad inter- 
est. Among them, the field induced transverse antifero- 
magnetism have been widely investigated in many ex- 
periemntal and theoretical studies.^""* If the magnetic 
field larger than the spin gap is applied to the spin 
gapped system, the single chain ground state becomes 
the Tomonaga-Luttinger liquid and the transverse anti- 
ferromagnetic order develops as soon as the weak inter- 
chain coupling is swtiched on. 

The disorder is another origin of order in the spin- 
gapped low dimensional quantum magnets. In the 
presence of disorder, the spins in the nonmagnetic ground 
state revive and induce so-called 'disorder induced order' 
even in the absence of the magnetic field. 

On the other hand, the possibility of disorder induced 
magnetization plateau is also predicted in a certain class 
of one-dimensional random quantum magnets. ^"^^^ This 
corresponds to the spin gap state induced by disorder 
and magnetic field. 

Therefore, the effect of disorder on the quantum mag- 
nets in the magnetic field is twofold. Namely, it en- 
hances the magnetic order by reviving the spins, while 
it suppresses the magnetic order by forming plateaus. 
In the present work, we investigate the competition be- 



tween these two controversial aspects of randomness in 
the quasi-one dimensional quantum spin systems and the 
resulting multiple reentrant phase transitions between 
the transverse antiferromagnetic phases and disordered 
plateau phases. The similar problem has been discussed 
in the diluted dimcr network system by Mikeska and 
coworkers^^ and in the coupled random dimer network 
by Nohadani and coworkers. 

This paper is organized as follows. In the next section, 
the model Hamiltonian is presented. The single chain 
magnetization curve is calculated in section 3. Within 
the interchain mean field approximation, we predict the 
multiple reentrant behavior with the increase of the mag- 
netic field in section 4. The calculation of the spin-spin 
correlation function is presented in §5. Even in the non- 
plateau state, where the infinitesimal interchain coupling 
induces the transverse ordering, the correlation function 
of a single chain turned out to be short ranged. Based on 
these observations, the mechanism of the antiferromag- 
netism away from the plateau region is explained in §5. 
The final section is devoted to summary and discussion. 

2. Model Hamiltonian 

As a candidate model in which the reentrant antifer- 
romagnetism is expected, we investigate the quasi-one- 
dimensional random dimerized S — 1/2 Heisenberg chain 
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Fig. 1. The magnetization curves of pure single chains with p=0 
and 1.0. The chain length is N = 240. 

whose Hamiltonian is given by, 

{N/2 N/2 1 

y^^JS2i-l,jS2i,j + ^ JijS2i,jS2i+l,j > 
i=l i=l J 

TV 

+ ^ ^ JintSijSij' , (1) 
»=1 <3J'> 

where Jj = Js with probabiUty p and Jj = Jw with prob- 
abiUty 1—p- The interchain exchange coupling is denoted 
by Jint • The spin operator Sij denotes the spin on the i- 
th site of the j-th chain. The summation X]<j j'> taken 
over all nearest neighbour pairs of chains. In the present 
work, we assume Js > J > Jw > 0. Similar model with 
ferromagnetic Jw has been discussed^^ related with the 
experimental materials.^ 

3. Single Chain Magnetization curve 

The ground state magnetization curves of the single 
chain with Ji„t = in (1) is calculated using the DMRG 
method. The magnetization per site M is defined by 

1 ^ 

where the summation is taken over all spins in a single 
chain and the chain index j is suppressed. The regular 
models with p = or p = 1 has magnetization plateaus 
at M = which corresponds to the spin gap and at the 
saturation magnetization M = Mg = 1/2. However, it 
has no plateau with intermediate values of magnetization 
as shown in Fig. 1. 

On the other hand, the magnetization curves for p ^ 
0, 1 consist of a sequence of plateaus. Between them, the 
magnetization increases almost continuously. The typical 
example is shown in Fig. 2 for Js = 2, Jw = 0.1 and 
J = 1 for various values of p. 

The main features of the magnetization curves can be 



Fig. 2. The magnetization curve of single chains for p=0.2, 0.4, 
0.6 and 0.8. The chain length is N = 120. The magnetization 
is measured for the middle 60 sites to reduce the boundary ef- 
fect. The average is taken over 64 samples. The error bars are 
shown only for selected points because otherwise the symbols are 
extremely dense. 

understood using the cluster picture similar to that de- 
scribed in ref.^° With the increase of magnetic field, fol- 
low three large plateaus which are numbered I, II and 
III in Fig. 2. Let us consider a cluster consisting of q 
successive Js-bonds and q — 1 J-bonds in between. This 
is called '(/-cluster' as in ref.^" The 2q spins in a cluster 
form a tightly bound singlet cluster. The two spins con- 
nected to both ends of this cluster by J-bonds are almost 
free but weakly coupled mediated by the quantum fluc- 
tuation within the strongly coupled cluster. Other spins 
form singlet dimers on the J-bonds if Jw « J- 

On the plateau I with magnetization M = (1 — p)Ms, 
the spins which do not belong to the g-clusters are all 
polarized. On the plateau II, the end spins separated 
by 1-clusters with a single Js-bond remain unpolarized. 
Similarly, on the plateau III, the end spins separated by 
2-clusters also remain unpolarized and so on. These in- 
terpretation are confirmed by comparing with the mag- 
netization process of a cluster consisting of a g-cluster 
and two additional end spins connected by J-bonds on 
both ends of the g-cluster. Lower plateaus due to the 
spins separated by longer g-chistcrs are not clearly iden- 
tified within the present scale. The low field part of the 
magnetization curve reflects the singularity of the low 
energy excitation spectrum as described in rcf. 10. 

Above the plateau I, the magnetization increases with 
series of plateaus and narrow continuous parts up to the 
saturation field. As p increases, the; width of the plateaus 
decrease and magnetization curve becomes almost con- 
tinuous. 

4. Effect of Interchain Exchange Interaction 

We treat the interchain coupling by the mean field ap- 
proximation^^ assuming the transverse antiferromagnetic 
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(3) 



order as, 

{-lym Jint < 

i-iyPjm Jint > 

For Jint > 0, we assume that the two dimensional lat- 
tice of the chains is bipartite. The quantity Pj is +1 if 
the site j belongs to one of the sublattice and —1 if it 
belongs to the other. We thus have the interchain mean 
field Hamiltonian H^^^ for each chain as 



N 



N 



H 



IMF 



J S2i-lS2i + JiS2iS2i+l 



1=1 



N 



(4) 



with Jfst = -z| Jint|m. 

In order to investigate the stability against transverse 
antiferromagnetic ordering, we fix the interchain mean 
field Hgt and calculate the finite field staggered suscepti- 
bility x(i?st) = m{Hst)/Hst which tends to the staggered 
susceptibility Xst as iJgt tends to 0. In general x7t^ gives 
the minimum interchain coupling A — 2;| Jint| which sta- 
bilizes the transverse ordering. If Xst diverges, the trans- 
verse ordering takes place for the infinitesimal interchain 
coupling within the interchain mean field approximation. 

Fig. 3 shows the magnetic field dependence of 
Xst\H,t) with Hst = 0.0005 for p = 0.2 and p = 0.8 
as representatives of the small p and large p cases. The 
magnetization curves in the absence of staggered field 
is also presented. The H-x~^{Hst) curve has multiple 
maxima, which clearly shows that multiple reentrant be- 
havior takes place for finite interchain coupling in the 
ground state. It should be noted that the H-x~^{Hst) 
curves are insensitive to the value of iJgt around these 
maxima. However, around the dips and minima, the val- 
ues of x{Hst) have significant /fst-dependence. 

For p = 0.8, x^^(iJst) remains significantly small 
above the main plateaus compared to the peak values 
on the main plateaus. Therefore we expect no disordered 
phase for moderate values of interchain coupling in this 
region where the magnetization curve appears almost 
continuous. Even away from the main plateaus, however, 
the magnetization curve shows a series of small plateaus 
and narrow continual parts between them. Correspond- 
ingly, x^^(iJst) tends to a finite value as Hst — » on 
these plateaus. The detailed features of such behavior is 
shown in Fig. 4 for 2.17 < H < 2.21 as a representative. 
On these small plateaus, x^^(-ffst) clearly tends to small 
finite values as shown in Fig. 4 around H = 2.173 and 
H = 2.208. On the other hand, in the true off-plateau 
state, Xst tends to zero suggesting the divergence of Xst- 
In this case, the transverse antiferromagnetic order is sta- 
bilized in the presence of infinitesimal interchain interac- 
tion. In Fig. 4, such behavior is observed at H = 2.1905. 
To investigate this behavior in more detail, we present 
the i?st-dependence of Xst in Fig. 5 at H = 2.1925, 2.1905 
and 2.188 using the data for N =480. Only 0.2% devia- 
tion from H = 2.1905 causes a clear upturn of x^'^(iJsi). 
Although the small size dependence is present, the data 
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Fig. 3. The i?-dcpondence of l/xst(Hst) with H^t = 0.0005 for 
(a) p = 0.2 and (b) p = 0.8. The magnetization curves are also 
shown for reference. The staggered magnetization m is measured 
for the middle 60 sites to reduce the boundary effect. The chain 
length is AT = 120 and average is taken over 512 samples. The 
error bars are shown only for selected points. 
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Fig. 4. The //-dependence of l/xstC-ffst) for 2.17 < H _ 
2.21. The values of the staggered field are Hst = 
0.002,0.0015,0.001,0.0005,0.00025 and 0.0001 from top to bot- 
tom. The chain length is iV = 480 and measurement is done for 
the middle 240 sites. The average is taJsen over 256 samples.The 
error bars are within the size of the symbols. 
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Fig. 5. The //st-dcpcndcncc of Xst{J?st)^^ around H = 2.1905. 
The big symbols arc for N = 480 and small ones are for N = 
240. The solid lines are power law extrapolation from Hst = 
0.0005, 0.001, 0.0015 and 0.002 for N = 480. The average is taken 
over 256 samples. 
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Fig. 6. The transverse correlation function | SJy | at = 
2.1905(o) where Xst(-ffst ^ 0) = and at H = 2.1735(») where 
Xst(^^st — * 0) = finite. The average is taJjen over 512 samples. 




Fig. 7. The log-log plot of the correlation function with p = 0.8 
(open circles) and p = l(filled squares) at H = 2.1905. 
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Fig. 8. Schematic finite temperature phase diagram on the H—T- 
plane. 



for N = 240 also shows the similar behavior as plotted in 
smaller symbols. Therefore we expect this is not the finite 
size eflFect but is an essential feature of the present system 
in the thermodynamic limit. Similar behavior is observed 
for other values of H above the main plateaus. Therefore 
we conclude that the antiferromagnetic order is stabi- 
lized by infinitesimal interchain coupling only within a 
narrow region where magnetization increases continously 
between the successive small pleteaus. 

5. Correlation Functions 

In order to get more insight into the nature of 
each state based on the properties of single chains, we 
have also investigated the spin-spin correlation function 
(^SfSj) as a function of |j — j|. Figure 6 shows | (5f 5J) 
for H = 2.1905(o) where Xst(-ffst) tend to zero. Even 
in this case, the spin-spin correlation function is short 
ranged. Actually, the behavior of the corelation func- 
tion is almost the same as that for H = 2.1735(») where 
Xst(-ffst) tends to a small but finite value. In Fig. 7, the 
log-log plot of the same correlation function is compared 
with that for the regular chain with p = 1 at H = 2.1905. 
It is clear that the rapid decrease of the correlation is dis- 
tinct from the power decay for the regular chain. 

This can be understood in the following way. In the off- 
plateau region, the continuum of the low energy excited 
states pile up on the ground state. In many of these ex- 
cited states, the spins which are not correlated in the 
ground state are correlated. The staggered transverse 
magnetic field mixes up these excited states and leads 
to the divergent staggered susceptibility. In this case, 
the long range transverse order can be stablized with 
infinitesimal interchain coupling even though the spin 
correlation is short ranged in the ground state. On the 
other hand, in the plateau state there exists no low en- 
ergy excited states which supports the long range order 
with small interchain coupling. 

6. Summary and Discussion 

The transverse magnetic ordering in the ground state 
of the random quantum Heisenberg chain is investigated 
using the density matrix renormalization group and the 
interchain moan field approximation. It is predicted that 
the multiple reentrant behavior takes place between the 
disordered plateau phases and transverse antiferromag- 
netic ordered phases. This is in contrast to the case 
of random dimer networks discussed by Mikeska and 
coworkers^^ and Nohadani and coworkers^^ for which the 
reentrant transisition takes place only once. 

It is also pointed out that even in the non-plateau 
regime the spin-spin correlation of the single chain is 
short ranged. Nevertheless, the long range order is estab- 
lished with infinitesimal interchain interaction with the 
help of the excited states which pile up near the ground 
state and are mixed up by the interchain staggered mean 
field. 

In this work we concentrated on the ground state 
phase transition. Nevertheless, the reentrant behavior 
should survive even at finite temperatures. The transi- 
tion temperature should be high between the plateau re- 
gion and low or zero on the plateau region as depicted 
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in Fig. 8 schematically. This behavior manifests itself as 
anomalous behavior even above the transition temper- 
ature. Therefore, if we increase the magnetic field with 
fixed temperature and passes near the ordered phase, the 
transverse spin fluctuation would be strongly enhanced. 

We expect the present type of reentrant behavior is 
universal in the random quantum spin systems in which 
the singlet dimer formation is randomly perturbed to 
produce the local almost free spins. In contrast to the 
random dimer network systems in which the dilution pro- 
duces the isolated spins, the free spins in the present sys- 
tem are produced by the random competition of two dif- 
ferent dimer interactions Js and J each of which prefers 
different dimer configuration. This is the origin of more 
complicated structure of the phase diagram. Thus we ex- 
pect the reentrant behavior of due to similar mechanism 
in a variety of systems. 

The computation in this work has been done using 
the facilities of the Supercomputer Center, Institute for 
Solid State Physics, University of Tokyo and the Informa- 
tion Processing Center, Saitama University. This work is 
supported by a Grant-in- Aid for Scientific Research from 
the Ministry of Education, Culture, Sports, Science and 
Technology, Japan. 
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